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ABSTRACT 

Aims. A model of heliosheath density and energy spectra of a-particles and He + ions carried by the solar wind is developed. 
Neutralization of heliosheath He ions, mainly by charge exchange (CX) with neutral interstellar H and He atoms, gives rise to 
~0.2 — 100 keV fluxes of energetic neutral He atoms (He ENA). Such fluxes, if observed, would give information about plasmas in 
the heliosheath and heliospheric tail. 

Methods. Helium ions crossing the termination shock (TS) constitute suprathermal (test) particles convected by (locally also diffusing 
through) hydrodynamically calculated background plasma flows (three versions of flows are employed). The He ions proceed from the 
TS towards heliopause (HP) and finally to the heliospheric tail (HT). Calculations of the evolution of a- and He + particle densities and 
energy spectra include binary interactions with background plasma and interstellar atoms (radiative and dielectronic recombinations, 
single and double CX, stripping, photoionization and impact ionizations), adiabatic heating (cooling) resulting from flow compression 
(rarefaction), and Coulomb scattering on background plasma. 

Results. Neutralization of suprathermal He ions leads to the emergence of He ENA fluxes with energy spectra modified by the 
Compton-Getting effect at emission and ENA loss during flight to the Sun. Energy-integrated He ENA intensities are in the range 
~0.05 50 cm~ 2 s sr -1 depending on spectra at the TS (assumed kappa-distributions), background plasma model, and look direc- 
tion. The tail/apex intensity ratio varies between ~1.8 and ~800 depending on model assumptions. Energy spectra are broad with 
maxima in the ~0.2 — 3 keV range depending on the look direction and model. 

Conclusions. Expected heliosheath He ENA fluxes may be measurable based on the capabilities of the IBEX spacecraft. Data could 
offer insight into the heliosheath structure and improve understanding of the post-TS solar wind plasmas. HT direction and extent 
could be assessed. 

Key words. Sun: heliosphere - particle emission - Plasmas - Atomic processes - Accelerations of particles - ISM: atoms 



1. Introduction 

Helium, as the second abundant ion species in the solar wind, 
should also be prominent in energetic neutral atom (ENA) fluxes 
from the heliosheath, resulting from the transcharge on neutral 
atom populations. Detection of these fluxes is contingent upon 
the energy of the He ENAs. It is known that He ions at in- 
terplanetary shocks do not equilibrate their downstream ther- 
mal energies with protons. SWICS experiment data obtained on 
board Ulysses s/c for 15 shocks at distances between 2.7 and 5.1 
AU from the Sun indicate that the mean ratio of a-to-protons 
down stream thermal velocities is 1.3 + 0.3 (Berdichevsky et al. 
1997). High thermal speeds of c-particles (and pickup He + ) at 



tion time for 1 keV/n a-particles (i.e. bulk pre-TS He) with he- 
liosheath thermal proto ns of number d ensity ~0.001 cm 3 and 



shocks accompanied by enhanced high-energy power-law tails 
were also seen by Ulysses SWICS and HISCALE experiments 
in a reverse quasi-perpendicular CIR shock and also at the in- 
bound pass of the quasi-parallel Jovian bow shock (Gloeckler 
|et al.|200"5j ). Instances of non-thermalization of He ions at Earth 
bowshock were reported and discussed by |Liu et"aL ( 2007) 1. 

All this suggests that the He ions in the heliosheath can be 
treated at injection, i.e. immediately upon crossing the TS, as 
thermally separate from the bulk of shocked solar wind plasma: 
the average post-shock energy should be close to the upstream 
energy in a solar frame, with tail distributions approximated by 
power laws in energy. The excess energy, compared to ideal gas 
dynamics, of the random motions of He ions in the heliosheath 
can probably persist for a long time. The Coulomb equilibra- 



2011 



amounts to 5 x 10 s, 



temperature ~70 000 K (Richardson ; 

which is longer than the heliosheath residence time ~10 9 s. Loss 
of energy through interaction with neutral H atoms in the he- 
liosheath is also slow. For total stopping p ower for 1 keV/n a- 
particles in H equal to 1120 MeVcm 2 g _1 (Berger et al. 2005 1 
the characteristic time for c-energy decay in neutral H density 
0.2 cirT 3 is ~2 x 10 11 s. All this suggests that He ions in the he- 
liosheath will constitute a relatively hot particle population, with 
typical energies ~0.5-l keV/n, so well above the thermal proton 
plasma measured by Voyager, and with nonthermal tails reach- 
ing several keV/n. Neutralization of He ions, mainly by charge 
exchange with neutral interstellar H and He, should therefore 
lead to the emergence of fluxes of relatively energetic He ENAs. 
The aim of the present paper is to estimate the magnitude of the- 
ses fluxes and to look for the conditions for detection. 



Our approach is in part similar to the one used by 
Grzedzielski et al. (2010) to study the fate of solar-wind heavy 
ions in the heliosheath. The He ions are treated as test particles 
undergoing various binary interactions (BI), with other particle 
populations constituting the heliosheath plasma. This allows us 
to calculate (1) how the He ions change their charge states and 
energies, (2) what the emissivity of He ENA is as a function 
of energy for each point in the heliosheath, and (3) how many 
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of these ENAs, having survived losses on the way, will reach 
Earth's vicinity from a particular direction. 

2. Physical model 

In our model we follow the time evolution of charge states and 
energy of He ions in the heliosheath as solar wind plasma flows 
from the TS, towards the heliopause (HP), and finally to the he- 
liospheric tail (HT). We treat a-particles and He + ions as test 
particles, carried by the general flow, which undergo BI with 
background electrons and protons, with solar ionizing photons, 
and with neutral H and He atoms coming from interstellar space. 
In addition we take He ions energy changes into account due 
to (adiabatic) compression/decompression of the background 
flow, as well as energy change (in fact decrease) resulting from 
Coulomb scattering on plasma background and possible effects 
of spatial diffusion. 

2.1. Evolution of helium ions in phase space 

To describe the behavior of the two He ion species, a-particles 
and He + ions, we calculate changes in the local velocity dis- 
tribution functions f° and / He resulting from displacement of 
the considered plasma parcel along its flow line determined 
by the hydrodynamic time-independent solution for background 
plasma. The background flow is assumed to be stationary in time 
and axially symmetric, depending on distance r from the Sun and 
angle 9 from the apex axis (cf. Sect. 2.3 1. The functions /" and 
/ He are assumed to be isotropic in velocity space; that is, they 
are functions of the scalar momentum p. Their dependence on r 
and 8 can be expressed as functions of the curvilinear coordinate 
s along the flow line or, equivalently, using 



ds 



v\dt, 



(1) 



as functions of time f, i.e., the flow history of the parcel of back- 
ground plasma. v sw describes the solar wind bulk velocity in the 
heliosheath. 

The changes in f" and f He along the flow line are thus de- 
termined by coupled equations of the type used to describe the 
transport of cosmic rays (Joki pii|1987 1 



as 



-I* Wfe - L c a\p - L aM + V • (k • V/ 1 *) , (2) 



Equations Q and ([3} do not contain Fokker-Planck type 
terms that would describe possible local stochastic acceleration. 
This acceleration is often invoked in the context of ACR popula- 
tions, though its relevance - in view of the VI and V2 data - is 
debatable (Florinski et al., 2011, "The global heliosphere during 
the recent solar minimum" talk at the Solar Minimum Workshop, 
Boulder CO, May 17 - 19, 201 1). Applying stochastic accelera- 
tion to the present context of ~ keV ions would require a sound 
understanding of the small-scale magnetohydrodynamic turbu- 
lence in the heliosheath, which is lacking at present. 

In numerical solution of Eqs. Q and ^ we calculate the 
evolution of a and He + spectra separately for each flow line. We 
set initial discretized spectra at the TS (cf. Sect. 2.5 1 and then 
calculate the evolution between adjacent points. Discretization 
is fixed at 500 bins between km s _1 and 5000 km s _1 with a 
10 km s width each. As the plasma parcel proceeds along its 
flow line, particles are shuffled between bins as required by the 
interactions. This scheme is repeated for all flow lines. 

To be able to calculate the righthand side of Eqs. Q and <|3j, 
one should know, besides the relevant cross sections, the back- 
ground solar plasma and neutral interstellar gas flows, e.g., the 
solar wind electron (n e ), proton (n p ) densities, the density dis- 
tribution of interstellar neutral H and He (nu, «He), as well as 
the corresponding bulk flow velocities of the heliosheath solar 
wind plasma (v sw ), interstellar gas (vh), and the effective relative 
velocities of particles at collisions (v re i) resulting from local par- 
ticle velocity distribution functions. Modeling of these functions 
is described in Sect. 12.31 



„ f He+ - 1 T7 „ 9 fHe + , r T I 

di _ 3 V sw 0(\np)J + "BI,mHe* _ ^BI,He + ^ff + 

— ^BI,He + — >He ~ L C ,vf +Xp ~ ^He+,H + V • \K • V/ He ) (3) 

The successive terms on the righthand side of Eq. Q describe 
changes in f a due to adiabatic compression/rarefaction of the 
background flow, gain (GBi,He + ^a) from BI conversion of He + 
into a, loss (L Bla ^ He +) from BI conversion of a into He + , and 
loss (Lfii.a^He) from BI conversion of a into He, loss (L c v „\ P ) due 

to Coulomb scattering on backgro^rotons corresponding to £ 3 Heljosheath background p/asma 



2.2. Binary interactions affecting the He ions 

We assume that the He ions (test particles) are immersed in a 
substratum constituted by background heliosheath protons and 
electrons, background neutral H and He atoms of interstellar ori- 
gin, and ionizing solar photons. The binary interactions include 
radiative and dielectronic recombinations, electron impact ion- 
izations, photoionizations, double and single charge exchanges 
(also to upper levels), and electron stripping. The scheme of tran- 
sitions between the charge states of He ions resulting from the 
interactions is shown in Fig.[T] in which the three levels describe 
He charge-states 0, +1, +2 and the arrows correspond to binary 
interactions denoted a, b, c, d, e, f, g, h, i, k, 1, m (Table [TJ. 
As we look for He ENAs with velocities > 100 km s _1 (energy 
> 207 eV), we take the He + pickup ions born in the supersonic 
solar wind into account, while we disregard those originating in 
the heliosheath, where the relative velocity between heliosheath 
flow and interstellar neutral He atoms is rather low. 

Interactions b, c, d, e, k, 1, m shuffle the ions between charge 
states 1 and 2. They preserve the total number of ionized He 
atoms. Interactions a, f, g, h, i convert He ions into He ENAs 
and are the source of presumed He ENA fluxes at Earth. Once a 
He ENA is born, it is assumed to be lost to the heliosheath He 
budget. Probability of reionization of a 1 keV He atom flying 
from 5000 AU in the tail to Earth is only 12%. 



Huba 



( 2002 1, and loss (L a , H ) 



energy loss rate v e as given by 
due to a interaction with neutral hydrogen (Berger et al. 2005). 
The last term on the righthand side describes possible effect of a- 
particles spatial diffusion, with corresponding tensorial diffusion 
coefficient k (cf. Sect. |2.4[ ). Equation Q is analogous to Eq. (|2]i 
with symbols a and He + interchanged. Details on BI are given 
in Sect.l2~2l 



The state of bulk solar wind plasma in the heliosheath is, at 
present, the subject of an intense debate brought about by un- 
expected results of the plasma experiment on board Voyager-2 
(V2) (Richard son et al.||2008p during and after crossing(s) of 
the TS in Aug./Sep. 2007. There is little doubt that the down- 
wind (post-TS) plasma is in a very different state from what was 
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Fig. 1. Scheme of transitions between charge 
states of He resulting from the considered bi- 
nary interaction with heliosheath background 
ion, electron, and neutral atom populations (cf. 
text). The heavy arrows denote transitions of 
greatest importance in most of heliosheath re- 
gions, the dashed arrow (crossed over) is the 
neglected production of low-energy He + ions. 
Emissivity (source funtion) of He ENA is de- 
termined by transitions a+f+g+h+i. 



Table 1. Binary interactions determining charge-state changes 
of He ions in the heliosheath 



Sign 


Reaction 


Reference 


a 


ct, He — » He, a (double cx) 


(1) A104 v 


b 


a radiative recombination with e 


(2) 


c 


a, H -> He + , p 


(3) z 


d 


a, He -> He\ He + 


(1) A98, B112 


c 


or, He + -> He + , a 


(1) A84 


f 


He + radiative recombination with e 


(4) 


g 


He + dielectronic recombination with e 


(4) 


h 


He + ,H-> He,p 


(1) A62 


i 


He + , He -> He, He + 


(1) A70 


k 


photoionization of He + 


(5) 


1 


electron impact ionization of He + 


(6) 


m 


stripping/cx ionization by p 


(1) A48, A54 



Notes, e - electron, p - proton, a - a-particle, cx - charge exchange 
w A 104 and similars symbols denote reactions as listed in (1) 
<z) (1) A88 for energy > 5 keV, excitation and/or emission involving 
upper levels includes (1) B90, B96, B102, B104 



Re ferences. (1) |Bamett| <fT990)l; (2) |Arnaud & Rothen flug] (1985); 



(3) Liu et al. (2 003); (4)|Aldrovandi & Pequignof[ p973) ; (5) |BochsIer| 
|et al.H2012) ; (6) |Janev et al-H1987> 



expected on the basis of standard Rankine-Hugoniot equations. 
The post-shock temperature of the majority (thermal) protons 
seems to be much lower (~70000 K) than expected in a single- 
fluid shock transition (~10 6 K), and the bulk flow velocity starts 
to decrease well ahead of the shock with a much smaller velocity 
jump at the shock itself. Also the electron temperature seems to 
be quite low, T e < 10 eV ( |Richardson|2008| l. 

In contrast, a relatively small fraction (~10-30%) of post- 
shock protons is endowed with energies of ~ 1 to perhaps several 



keV. This nonthermal proton population, resulting presumably 
from ionized interstellar H atoms picked up by the supersonic 
solar wind, is thought to contain the bulk of total energy den- 
sity (pressure) of the post-shock plasma. V2 data indicate such 
a situation seems to prevail deep into the post-shock plasma 
(Richardson 2011) with little spatial change observed in ther- 
mal proton density and temperature. Such conditions may ex- 
ist in the frontal lobes and perhaps also in near tail of the he- 
liosheath. Another important fact is a much faster than antici- 
pated decline in heliosheath plasma bulk velocity derived from 
Voyager- 1 (VI) data as the spacecraft receded from the TS in 
years 2005-2011 ( |Krimigis et al. 201 1) . In particular the si- 
multaneous drastic decrease on Aug. 25, 2012 in the fluxes of 
ions of energies > 0.5 MeV/n observed by the LET telescopes 
of the cosmic ray subsystem instrument on board VI and the 
accompanying increase in the magnetic field to about 0.41 nT 



(ftp://lepvgr.gsfc.nasa.gov/pub/voyager/) seem to indicate that 
VI might have crossed the heliopause around this date. 

These results imply that the TS-HP stretch along VI trajec- 
tory is only ~ 121.5 - 94 =~ 27.5 AU long; that is, the he- 
liosheath is much narrower than ~60 100 AU obtained in 

gasdynamical modeling based on data available prior to VI and 
V2 TS-crossings (e.g., in the Izmode nov & Alexashov| ([2003 ) 
model the distances Sun-TS and Sun-HP for ecliptic latitudes 
corresponding to the V2 trajectory are 110 AU and 208 AU, re- 
spectively). A narrow heliosheath of only 25 ± 8 AU thickness in 
the upwind region also comes out from the analysis of H ENA 
fluxes obse rved by the IBEX SOHO/HSTOF, and Cassini/INCA 
spacecraft ( Hsieh et al.|2010 1. 

These new findings are at present not properly integrated and 
understood within a coherent physical picture. Despite the suc- 
cess of a realistic and time-varying description of the TS cross- 
ing positions by both Voyager spacecraft in a recent model by 
Washimi et al. (201 1|), the single-fluid MHD calculations em- 



ployed therein are unable to render the nonthermal aspects of the 
particle distribution functions in the heliosheath that may be im- 
portant for He ions physics. Therefore to numerically describe 
the background plasma flow conditions we employ three time- 
independent, axisymmetric heliosheath models (denoted in the 
following: hydrodynamic, Parker, and ad hoc), corresponding to 
three simple variants of assumed heliosheath plasma. 

The hydrodynamic model is the model developed by 

and 



Izmodenov & Alexashov 
Grzedzielski et al.| ( POLO) ). 



(2003 



The 



used previously in 
'arker and ad-hoc models, 



though not internally coherent in terms of physics, are 'tailored' 
in such a way as to approximately render the V2 suggested spa- 
tial distribution of variables decisive for the behavior of He ions: 
that is, thermal proton and electron densities (n p , « e ), tempera- 
tures (T p , T e ), as well as the main traits of the nonthermal plasma 
components. We also try to approximately render the common 
heliosheath bulk flow velocity (v sw ) along the V2 trajectory. For 
each variant, integration of Eqs. (|2]( and ([3]) is performed. In this 
way we follow the time evolution of the charge states, spatial 
distribution, and energy spectra of He ions in each fluid element 
carried by the background flow. This allows the local He ENA 
production rates to be calculated and - after accounting for He 
ENA energy losses and reionization on the way - we can con- 
struct expected He ENA spectra at Earth. The point is to test 
how sensitive the predicted He ENA fluxes are to the assumed 
widely discordant variants of plasma. We believe that if impor- 
tant common features of predicted He ENA fluxes appear in all 
considered cases, credence could be given to the results despite 
the partial inadequacy of the physical modeling employed. 

In the hydrodynamic model, the background flow of solar 
plasma and neutral hydrogen atoms in the supersonic solar wind, 
inner heliosheath, and distant heliospheric tail is described as 
single-fluid, non-magnetic, gas-dynamical flow of heliospheric 
proton-electron plasma coupled by mass, momentum, and en- 
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ergy exchange with the neutral interstellar hydrogen atoms cal- 
culated kinetically (Monte-Carlo approach). The Sun as a source 
of solar wind and ionizing photons is assumed to be spherically 
symmetric, with the wind speed of 450 km s _1 , Mach num- 
ber 10, and n p = 7 cirT 3 at Earth orbit. At infinity, a uniform 
interstellar flow of 25 km s _1 with neutral hydrogen density 
«h,lism = 0.2 cm" 3 , proton density holism = 0.07 cirT 3 , and 
temperature 6000 K was assumed. 

Interstellar neutral He atoms are represent ed by a uniform 
substr atum with density «He = 0.015 at. cm 3 (Gloeckler et al. 
2004) flowing with velocity of 25 km s _1 . Because of axial sym- 
metry, all variables depend on the radial distance r from the Sun 
and angle 8 from the apex direction. The background flow is 
found in form of n p , n e , T p (= T e ) given as functions of dis- 
tance s along individual flow lines. There are 180 flow lines, 
each starting at Earth's orbit. The flow lines are identified by the 
initial (1 AU) value of the angle 8. In this solution the TS and 
HP are, correspondingly, 102 and 177 AU distant from the Sun 
along the apex axis. 

The shocked solar wind plasma is very hot immediately be- 
hind the TS (T p = T e ~10 6 K) and then cools down to about 
20000 K in the distant heliotail. The background plasma veloc- 
ity distributions are always local Maxwellians. 

The Parker model is the classical subsonic solution for a 
point source of incompressible fluid (stellar wind) in a uniform, 
incompressible external flow ( |Parker| fl961 ). To be applicable, 
the model requires the TS radius much less than the distance 
Sun-HP. We assume that the solar wind mass source and the 
interstellar flow are the same as in the hydrodynamical model. 
The background velocities and shape of the flow lines are de- 
termined by the Parker analytical solution. Typical heliosheath 
velocities are in the range 120-250 km s _1 at ~100 AU and 
tend to 170 km s in the tail. We crudely approximate ob- 
served heliosheath plasma conditions by assuming constant he- 
liosheath proton density to be constituted by two proton popu- 
lations, the thermal one with n p t h = 0.0015 cirT 3 and another 
nonthermal one with M pn th = 0.0005 cm 3 . The proton ener- 
gies are described by a Maxwellian with 70 000 K temperature 
and a monoenergetic population of 1.1 keV per proton, respec- 
tively. Electron temperature is set to T e = 3 eV. To conform 
to a smaller HP distance as presently envisaged (Krimigis et al. 
201 1\ the Sun-HP stretch along the apex axis is set to 101 AU. 
The corresponding Sun-HP distance at the V2 ecliptic latitude 
is 83.7 + 27.5 = 111.2 AU. Comparison of proton density, bulk 
radial, and tangentional velocities measured in the heliosheath 
by V2 with the Parker model solution is shown in Fig. [2] 

The ad hoc model was developed for this research. We in- 
clude in it both the new (smaller) heliospheric scale and some of 
the new physical aspects of heliosheath plasmas. The solar wind 
mass source and the interstellar flow are again the same as in the 
hydrodynamic model. However, the background flow in the he- 
liosheath is solved anew. We use the approximation that the flow 
along each of the flow lines is one-dimensional in a known chan- 
nel of varying cross section A(s). The geometry of the channels 
is determined by the flow lines of the hydrodynamic model but 
linearly rescaled so as to place the termination shock (with shape 
conserved) at the observed distance of the V2 crossing, that is, at 
83.7 AU from the Sun. In this rescaling, the heliopause is put at 
111.2 AU (at V2 ecliptic latitude). The flow along the reshaped 
flow lines is found by integrating equations of conservation of 
mass, momentum, and energy with the appropriate source terms 
on the righthand side. 



In a channel we have the mass conservation equation 



1 d 

Ads 



[^PpVsw] = Q , 



(4) 



where p p is total (thermal + nonthermal) mass density of pro- 
tons and v sw bulk plasma velocity in the heliosheath, and Q is 
the mass-loading term. We neglect the contribution of solar wind 
He ions and neutral interstellar He atoms to mass loading com- 
pared to the H contribution, and we also neglect electron impact 
ionization of H atoms because of low heliosheath electron tem- 
perature, r e = 3 eV. Mass loading results only from proton - 
H-atom charge exchange and net Q = 0. 

Momentum loading is nonvanishing, and we have in the 
equation of motion the frictional force (acting in direction s) re- 
sulting from proton - H-atom charge exchange: 



1 d 

Ad~s 



[ A Pp V sw] 



ds 



Pnth + 0"(v eff )VeffnHPp (v H cos 8 - v sw ) , (5) 



where o~ denotes the charge exchange cross section, and v e ff is 
the relative velocity between particles. In the energy equation 
we treat as small the kinetic energy of the bulk heliosheath flow, 
as well as the pressure of the thermal protons (kept at r Pitn = 
70 000 K), and retain solely the terms proportional to the high 
pressure (energy) p n ±, of the mono-energetic nonthermal proton 
population. We assume this energy (1) is convected with flow 
velocity v sw and (2) decays by charge exchange with interstellar 
H atoms on a time t cx depending on particle energy. Then the 
equation takes the form 



1 d 
Ad~s 



[Av sw yp n th] = (y- l)v s 



sw ! Pnth 

ds 



Pnth 



(6) 



where y is the usual adiabatic exponent. In numerical calcu- 
lations we assume the nonthermal population contains 25% 
of mass and average energy per proton at the TS is 1.1 keV 
(Giacalone & Decker 2010). This means that nonthermal pro- 
tons contain 77% of the total energy associated with particle 
random mot ions. The bulk velocity on the TS downwind side is 
150 km s _1 ( |Richardson et al.|2008[ ) and the TS downwind den- 
sity dependence on angle 8 is like in the hydrodynamic model 
but rescaled to the new TS position. The fit of the ad hoc model 
solution to V2 measurements in the heliosheath for proton den- 
sity, radial, and tangential velocities is shown in Fig. [3] 

The density in the ad hoc model seems to be too high by a 
factor ~2. A number of causes could be responsible, among them 
time-dependent effects in the solar wind flow (our modeling al- 
ways uses an 'average' solar wind model) and/or heliosheath 
asymmetry resulting from a skewed interstellar magnetic field. 

2.4. Adiabatic heating/cooling versus spatial diffusion 

Changes in background plasma density along the flow line in- 
duce He ion energy changes that can be treated as adiabatic if 
diffusion is slow enough. We made an estimate of the possible 
role of diffusion in the case of the hydrodynamic model. To find 
regions in the heliosheath in which diffusion of a-particles is 
negligible we use a simplified version of Eq. |2]) treated as the 
transport equation for pressure p a of the a-particle cosmic ray 
gas (Drury & Voelk 1981 1. In this approximation Eq. d2]i be- 



Pa 

dt \y- 1 



Pa 



7 Pa dp 
y - 1 p dt 



(7) 
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Fig. 2. Comparison of heliosheath proton density (left panel), radial velocity (middle panel), and tangential velocity (right panel) in 
the Parker model with corresponding in situ measurements ( Richards on|20 1 1 ) along the orbit of Voyager 2. 
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Fig. 3. Comparison of heliosheath proton density (left panel), radial velocity (middle panel), and tangential velocity (right panel) in 
the ad hoc model with corresponding in situ measurements (Richards on|201 \\ along the orbit of Voyager 2. 



in which p is the background plasma density. We assume dif- 
fusion is negligible as long as the first term on the righthand 
side (effect of diffusion) is much less than the second term (adi- 
abatic heating/cooling). The heliosheath diffusion coefficient for 
low-energy a-particles is in fact unknown. It depends most prob- 
ably on heliosheath turbulence. Models for the superposition of 



slab and 2-D turbulence were developed for the solar wind (Zank 
et al. 2004), but it is unclear how the assumed (2:8) energy par- 



tition between the two modes of turbulence corresponds to real 
heliosheath conditions. Therefore we take simply two limiting 
formulae as crude guesses: (i) Bohm diffusion, Kb, and (ii) phe- 
nomenological diffusion by |le Roux et al. (1996]) originally de- 
veloped for the ACR ions, kacr- We extrapolate them to our a- 
particle energies. In both formulae we assume k is scalar. 

Because diffusion coefficients (i) and (ii) differ, for the range 
of energies discussed (~0.2 - 20 keV, cf. Sect. [3]) by 4.5 to 5.5 
orders of magnitudes, we also take for comparison (iii) an inter- 
mediary diffusion coefficient equal to (kbkacr) 1 ' 2 - In calculating 
the values of the first term on the righthand side of Eq. (j7]) along 
the flow line, one has to use the local value of the heliosheath 
magnetic field B n sh- As the detailed structure of B n sh is unknown, 
we calculate B ns h along the flow line assuming the field starts on 
the downwind side of the TS with a value of 0. 1 nT as measured 
by |Burlaga et al. ([2005 2008 ) (variation in the magnetic field 
along the TS surface was described by standard dependence on 
angle/distance). Further evolution of the field then followed from 
background plasma density variations. In this we assumed that 
the magnetic field is frozen to background plasma, its direction 
is randomly oriented and plasma compression/decompression is 
isotropic, i.e. fihsh K (background plasma density) 2 ^ 3 . In this way 
one can construct maps of the ratio of first-to-second terms on 
the righthand side of Eq. (j7]( in the heliosheath. The maps are 
shown for values of the said ratio equal to 1/100 and 1/10 in 



Fig. [4] They correspond to k — and k = (kbaacr) 1 ^ 2 an d to 
a-particle energies at inject ion, peaking at energy correspond- 
ing to 170 km s _1 (cf. Sect. 2.5 1. Figure ^suggests that neglect 
of diffusion is justified for both diffusion versions (i) and (iii). 
For version (ii) diffusion is a paramount effect, but this case is 
hardly realistic for the low energies discussed here. Diffusion 
effects obtained for He + ions are similar. 

In the Parker model diffusion does not exist because of as- 
sumed uniform background plasma density. In the ad hoc model, 
the density distribution is quite flat (cf. Fig. [3] left panel), so dif- 
fusion is also negligible. In solving Eqs. (|2]i and ^ we therefore 
always consistently retain the adiabatic cooling/heating term 
(first on the righthand side). Energy variations due to these terms 
are moderate; for instance, energy increase along the flow line 
starting at 9 — 5° at the TS from apex, does not exceed a factor 
1.5 at its maximum. 



2.5. Density distributions and energy spectra of a-particles 
and He + ions in the heliosheath 

The a-particles (n a ) and He + ions («He + ) density distributions 
and energy spectra are calculated for the three background flow 
models (hydrodynamic, Parker, ad hoc, Sect. \2.3\ by integrat- 
ing Eqs. (|2]i and ([3]) along 180 flow lines. The approximations 
and procedures we adopted are described in Sects. 2.2 2.3 and 
|2.4| Consistent treatment applies effectively to the interval (100, 
2000) km s corresponding to energy range 0.207 - 82.9 keV 
for He ions. The initial values of n a and « H e + are stated at the TS 
for the hydrodynamic and ad hoc models, assuming that the total 
He content of the solar wind constitutes 5% of the local proton 
plasma (by number). For the Parker model they are stated at the 
distance fitted to TS crossing by V2. 
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Fig. 4. Heliosheath map of regions (light-gray 
areas) where in the hydrodynamic model the 
ratio of the lst-to-2nd terms on the righthand 
side of. Eq. l|7j is < 1/100 (lower panels) and 
< 1/10 (upper panels); that is, the role of dif- 
fusion is small. Left column is for Bohm dif- 
fusion, case (i); right column for (^b^acr)'' 2 , 
case(iii). 



Two variants of injected (initial) a-particle energy distribu- 
tions at the TS were taken: one peaking at energy 0.6 keV, which 
corresponds to velocity 170 km s _! (in the following we label 
this variant "170 km s '), the other peaking at energy 2.1 keV, 
which corresponds to 320 km s _1 (we label this variant "320 
km s _1 "). These energies were chosen because the first one could 
correspond to the case when excess post-TS energy of an a- 
particle is coming from the bulk velocity jump at the TS, and the 
second to the case when excess energy is approximately equal 
to total bulk kinetic energy at the shock (Richar dson et al.|2 008 ) 



(cf. Sects. [T]and 2.3 1. In both variants the initial velocity dis- 



tribution shape is the same: a kappa-distribution with k = 2.5. 
This value was chosen because it corresponds to the borderline 
between the near-equilibrium region applicable to heavy solar 
wind ions and the far-equilibrium region corresponding to the 
inner heliosheath (Livadiotis & McComas 201 1 , cf. their Fig. 2). 

The initial population of He + at the TS is made up of (a) 
solar wind a-particles singly-decharged during their flight to 
the TS and (b) pickup He + produced by neutral interstellar He 
ionization in the supersonic solar wind region. Based on solar 
wind as in the hydrodynamic model, we estimated that group (a) 
amounts to 0.0005 (by number) of the proton content. For group 
(b) the corresponding number is 0.002 (Rucinski et al. 2003). 



Initial spectra for He + group (a) at the TS are the same as those 
for a-particles. For He + group (b) we consistently use only one 
initial spectrum: a kappa(=2.5)-distribution that peaks at 4 keV. 

Energetic He ion heliosheath density distributions and en- 
ergy spectra obtained from integration of Eqs. (|2]) and Q are 
shown in Fig. [5] for variants "170 km s -1 " and "320 km s -1 ". 
The lefthand column of panels corresponds to "170 km s ', the 
righthand one to "320 km s _1 ". Each row corresponds to one 
of the three models employed (top to bottom: hydrodynamic, ad 
hoc, Parker). Upper (lower) half of each heliosphere map refers 
to a-particles (He + ions). To give some feeling of the evolution 
of energy spectra, above and below each map we show four pan- 
els with ion energy spectra corresponding to spatial positions 
labeled I, II, III, IV. As indicated, all these points lie on the same 
flow line starting at the TS at a point 90° away from the apex as 
seen from the Sun. 

Figure [5] illustrates some of the basic physics of He-ion be- 
havior. First densities per energy interval are much lower for He + 
than for a-particles. This reflects the differences in the injection 
at the TS. However, second, while the a-particles shift to lower 
energies when we go from point I to IV with monotonically de- 
creasing spectra preserved, the He + ions often develop maxima 
at mid-energy range, even though they decline at high energy. 

This has two main causes. One is that, on the whole, ener- 
getic He ions lose energy on average because of Coulomb scat- 
tering on the plasma background. Though adiabatic compres- 
sion may overturn this process, this happens only locally in the 



upwind heliosheath. In the end, both Coulomb scattering and 
general decompression towards the tail have to take over. The 
other cause is related to the circumstance that conversion of as 
into He + is, in practice, a one-way process. This comes from the 
high ionization potential of He + . Therefore while locally heated 
as may later resupply the mid energies of the He + population, 
no similar process operates in the opposite direction. 

Insight into the situation in the upwind heliosheath for the 
hydrodynamic model can be gathered from Fig. [6] which shows 
expanded maps of the frontal part of heliosphere in a similar for- 
mat to Fig. [5] The points I to IV are now shown along a flow 
line starting at a TS point 15° away from the apex. The en- 
ergy maxima of the He + ions are more pronounced here than 
in Fig. [5] This is due to more effective heating by compression 
in the frontal heliosheath. 

3. Expected He ENA fluxes from the heliosheath 

The local emissivity (source function) 7ena (cm _3 s _1 keV _1 ) 
of He ENA is determined by the product of reagents' densi- 
ties and reaction rates (i.e. relative vel ocity times cross sec- 



tion, cm s ) for BI a, f, g, h, i (Sect. 2.2 1. All these quan- 
tities can be found as functions of position in the three he- 
liosheath models (hydrodynamic, Parker, ad hoc) by procedures 
described in Sects. |TT]JT2j |0 |Z4j and [23] The He ENA in- 
tensity /ena(£, 6) (cirT^'sr^'keV -1 ) of particles of energy E 
coming from the direction at angle 6 from the apex is expressed 
as a line-of-sight (LOS) integral of the source function from the 
TS to HP 



/ena(£,0) = 



r 



jENA(E,r,B) 
4k 



F CQ (E,0,v m )X 

extO') 



exp 



•Jrrs 



dr' 



vena 



dr. (8) 



Formula |8]l contains corrections for re-ionization losses in 
the heliosheath and in the supersonic solar wind, as well as a 
correction (F cg) for the Compton-Getting effect. The losses in 
the heliosheath were calculated for BI as listed in Table UJ In 
Eq. ([H} they are described by the effective depth T ext for extinc- 
tion of He ENA with velocity vena- The main contributions to 
heliosheath losses come from reactions a, i, and s. In the hy- 
drodynamic model, He ENA reionization by electron impact is 
also of some importance in view of the high electron temper- 
ature. The ionization losses in the supersonic solar wind (L sw ) 
are mainly due to reaction p (photoionization) in Table [2] with 
a small contribution from reaction q inside the orbit of Mars 
( |Bzowski et al.|2012] >. On the whole, losses inside the supersonic 
solar wind are much less important than those in the heliosheath. 

The value of Fqg depends on He ENA energy E in the ob- 
server's frame and the velocity of the source carried by the back- 
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Fig. 5. Calculated density distributions and energy spectra of heliosheath energetic He ions for the hydrodynamic, Parker, and ad 
hoc models. Upper/middle/lower rows correspond to hydrodynamic/ad hoc/Parker models. Left column of the panels corresponds to 
the peak velocity of the initial kappa-distributions of a-particles injected at TS equal to 170 km s , right column, correspondingly, 
to 320 km s _1 . He + injection as described in the text. In a given row, the upper (lower) half of each heliosphere map represents 
density distribution of a-particles (He + ions). Distance scale is in AU. Density scales (equal for c-particles and He + ions only in the 
two upper rows) are indicated on vertical strips on the right side of each half-heliosphere map. Small panels placed above (below) 
maps illustrate the evolution of energy spectra of or-particles (He + ions) as the background plasma parcel is carried over points I, 
II, III, IV along the flow line. The line starting at the TS point 90° away from the apex is shown. Horizontal (vertical) axis in small 
panels describes ion energy in keV (ion density in cm^keV -1 ). 
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Fig. 6. Calculated heliosheath energetic He-ion density distributions and energy spectra for the frontal part of heliosphere (in the 
hydrodynamic model). Injection data in left panel correspond to "170 km s ', right panel to "320 km s '. Other details similar to 
Fig. [5] except that the flow line shown starts at a TS point 15° away from the apex. 



Table 2. Binary interactions determining He ENA losses in the 
heliosheath 



Sign Reaction Reference 



a : 


a, He — > He, a (double cx) 


(1) A104 


d : 


a, He -> He\ He + 


(1) A98 




He + , He -> He, He + 


(1) A70 


P 


He photoionization 


(2) 


q 


He electron impact ionization 


(3) 


r 


He, p — > He + , e, p 


(1) D24 


s 


He, p -> He + , H 


(1) A32 


t 


He, H — > He + , H 


(1) A10 



Notes, e - electron, p - proton, H - hydrogen atom, a - a-particle, cx - 
charge exchange; 
w same as in Table[T] 

References. (l) |Bamett| < [T990| >; (2) [Cummings et al.|(2002) ; (3) |Janev| 
[eTail ( fT987t 



model. The fluxes from upwind are smaller than those from the 
tail, but the difference is marginal for the hydrodynamic model. 
The decrease is only by a factor ~1.8. This is an effect of the 
He ion decharging relatively close to the Sun in the hydrody- 
namic model compared to other models. On the other hand, the 
tail/upwind contrast is significant for the Parker and ad hoc mod- 
els (ratios of ~60 and ~800, respectively). Detection from the tail 
should therefore be the easiest, and the contrast tail/upwind may 
tell something about background plasma conditions. 

Details on the dependence of He ENA energy spectra on 
model and angle 9 are given in Fig.[8]in which rows from top to 
bottom correspond to upwind (9 = 10°), crosswind (9 = 90°), 
and tail (9 = 170°) directions, for two variants of injection, 
"170 km s _1 " (left panels) and "320 km s _1 " (right panels). In 
all cases the spectra are broad in energy with maxima at ^ 1 keV 
for "170 km s _1 " and between ~0.5 to ~5 keV for "320 km s . 
Most energetic tail spectra are obtained in the ad hoc model. On 

the whole, the spectral region ~0.5 5 keV seems to be the 

most promising for detecting He ENA. 



ground relative to the Sun (v sw ) for a given model. Fcg is a func- 
tion of angular distance 9 of the LOS from apex. It is calculated 
for an observer at rest, at 1 AU from the Sun in the direction of 
incoming ENA, with solar gravity neglected. Fcg is particularly 
significant in reducing energy of He ENA coming from the up- 
wind heliosheath in the ad hoc model, in which background flow 
recedes quickly from the Sun. 

The dependence of total He ENA intensity (integrated over 
energy range 0.2 - 50 keV and with losses included) on angle 9 
from the apex is shown in Fig.|7]for the hydrodynamic, Parker, 
and ad hoc models for two assumed variants of ff-population 
injected at the TS ("170 km s _1 " and "320 km s _1 ") and He + 



injection as in Sect. 2.5 

Results presented in Fig. [7] indicate that expected energy- 
integrated intensities of He ENA seem to depend more on the 
background model than to be sensitive to the energy of He ions 
injected to the inner heliosheath plasma at the TS. However, 
in all three models the intensity from the tail is greatest and 
ranges from a few ENA (cm _2 s _1 sr _1 ) for the hydrodynamic and 
Parker models to a few tens of ENA (cm -2 s sr _1 ) for the ad hoc 



4. Final discussion and conclusions 

Calculating the heliosheath production of ENAs out of heavy 
solar ions that capture electrons from neutral interstellar atoms 
is basically a straightforward exercise, if ion densities are esti- 
mated, the flux of interstellar neutrals is measured, and particle 
collision velocities can be calculated based on simple assump- 
tions. This would apply if the state of background plasma were 
determined by hydrodynamics and if heavy ions (treated as test 
particles) were heavy enough to be weakly coupled to the back- 
ground and therefore were able to preserve energies they had at 
the TS. Such an approach has been used in the past (Grzedzielski 
|et al.|2010] > to follow de-ionization of multicharged coronal ions 
of C, O, N, Mg, Si, and S carried by the solar wind. 

In the present paper that deals with light ions such as a- 
particles and He + , we have been faced with the situation that, 
on one hand, the interaction with the background plasma can be 
more important and, on the other, the evidently non-single fluid 
state of the background (as established by Voyager) affects var- 
ious binary interactions to degrees that have not been assessed 
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Fig. 7. Total He ENA intensities (integrated over energy range 0.2 - 50 keV and with losses included) as a function of angle 8 from 
the apex for the hydrodynamic, Parker, and ad h oc models, and two variants of g-par ticle injection at the TS: "170 km s ' - left 
panel and "320 km s - right panel (Sect. 2.5 1. He + injection is as described in Sect. 2.5 



Table 3. Expected energy-integrated He ENA intensities 



Model 


hydrodynamic 


Parker 


ad hoc 




(tail - apex) 


(tail - apex) 


(tail - apex) 


Case "170 km s" 1 " 








intensity 


5-3 


5-0.1 


38-0.05 


(cm"~ 2 s -1 sr ) 








max. in spectrum 


0.3-0.7 


0.2 - 0.4 


1-0.2 


at energy (keV) 








flux (cm~ 2 s _1 ) 


25-21 


5-1.5 


95-5 


Case "320 km s _1 " 








intensity 


11-7 


15-0.2 


47-0.1 


(cm"" 2 s -1 sr ) 








max. in spectrum 


0.7-5 


1.3-1.5 


3-0.3 


at energy (keV) 








flux (cm~ 2 s _1 ) 


54-46 


16-4 


150-14 



reliably. This can make the outcome very sensitive to details of 
heliosheath modeling. 

In view of the present lack of understanding what goes on in 
heliosheath plasmas, we tried to explore the range of expected 
He ENA fluxes for some simple models covering a wide gamut 
of possibilities. If a feature has emerged in all models, one may 
presume it has some relevance. 



As explained in Sect. 2.3 we employed three simple descrip- 
tions of heliosheath plasma background: hydrodynamic model 
by Izmodenov & Alexashov (2003); classical Parker model; ad 
hoc model attempting to partially render the V 1 & V2 measured 
heliosheath data. These models cover a fairly wide range of pos- 
sibilities. For all three models we followed the evolution of a- 
particles and He + ions (isotropic) momentum distribution func- 
tions in phase space (Sects. |2~H|2.2||2~4] ), calculated heliosheath 
density distributions and energy spectra (Sect. 2.5 1, and esti- 
mated He ENA energy fluxes expected at Earth (Sect. |3j. The 
last estimates are essentially the observables given in Table [3] 
The data are organized following assumed peak injection energy 
of He-ions at the TS (cases "170 km s _1 " and "320 km s _1 ", cf. 
Sect.|23l 



IBEX is the first mission dedicated to the astronomy of en- 
ergetic neutral atoms (jMcComas 2009). It features two single- 
pixel, neutral atom cameras: IBEX- Hi ( jFunsten et al.|2009 1 and 
IBEX-Lo ( jFuselier et al.|2009[ >. Basically only IBEX-Lo has the 
capability of discerning species of the incoming atoms, includ- 
ing (indirectly) neutral He (Mob ius et"aLl|2009|l, and IBE X-Hi 
was designed to observe neutral H. |Allegrini et al.| ([2008 ) sug- 
gest an ingeneous technique by which some other species, in- 
cluding He, might also be registered owing to a special treatment 
of data from the IBEX-Hi anti-coincidence system. 

Results of this paper suggest that detecting maximum pre- 
dicted fluxes might be within the reach of IBEX instrumentation. 
The full possibilities of He ENAs measurements have not as yet 
been published, but the possibility of detection depends critically 
on the ratio of He ENA to H ENA fluxes ( |Allegrini et al.|2008). 
Based on published detection efficie ncies (|Allegrini et al.[2 008 



Fig. 2) and H ENA fluxes given by |McComas et aL ( 2010| > we 
preliminarily estimate that He ENA intensities predicted by our 
modeling could be measured from the heliospheric tail direc- 
tion alone. For the ad hoc model predictions, the measurements 
should be possible in the IBEX-Hi energy channels 1.1, 2.7, 
and 4.3 keV, both for the "170 km s" 1 " and "320 km s" 1 " vari- 
ants. For the hydrodynamic model predictions, the measurement 
might be marginally feasible at 2.7 keV for the "320 km s " 
variant. In all considered models, maximum He ENA flux should 
come from the tail, though the effect may be only marginal for 
the hydrodynamic model (cf. fluxes given in Table[3]). If an angu- 
lar tail-upwind decrease in the ENA signal could be measured, 
this might help to better pinpoint the essential physics of inter- 
actions. Maximum signal from the tail should also be related to 
maximum tail extension, which could help determine the direc- 
tion of heliosphere wake, independently of the IBEX hydrogen 
data as in McComas et al. (2012), with obvious inferences con- 
cerning the interstellar medium state. 

Since the expected He ENA intensity at 1 AU from the Sun 
strongly depends on the conditions in the inner heliosheath, He 
ENA detection, if accomplished, may become an important asset 
in diagnosing the physical state of the innner heliosheath plasma 
with embedded He ions, as a complement to diagnosing via H 
ENA currently performed by IBEX. 
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Fig. 8. Expected energy spectra of He ENA (cm _2 s _1 sr _1 keV _1 ) from the heliosheath. Left/right column - TS injection variants 
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